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Abstract 



The late time behavior of waves propagating on a general curved spacetime 

is studied. The late time tail is not necessarily an inverse power of time. Our 

work extends, places in context, and provides understanding for the known 

results for the Schwarzschild spacetime. Analytic and numerical results are 

in excellent agreement. 
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Introduction. Waves propagating on curved spacetime develop "tails" . A pulse of grav- 
itational waves (or other massless fields) travels not only along the light cone, but also 
spreads out behind it, and slowly dies off in "tails" jfflJ. This tail phenomenon is fascinat- 
ing theoretically, and has also been found in post-Newtonian calculations to have possibly 
observable secular effects on the phase of the orbit of inspiralling binary systems J7|. 

For asymptotically late times this tail often has a particularly simple behavior, namely, it 
decays as an inverse power in t. Detailed analyses have been carried out for the Schwarzschild 
geometry, using both analytic p|J|J^,|] and numerical techniques ||f|]. These works are based 
on the Regge- Wheeler perturbation formulation, in which the propagation of linearized grav- 
itational, electromagnetic and scalar waves is described by the Klein-Gordon (KG) equation 

D<P(x, t) = [d 2 t -d 2 x + V{x)] (f)(x, t) = 0, (f ) 

with V(x) describing the scattering of <fi by the background geometry. The late time tail 
has been explained in two different ways: in terms of a branch cut in the Green's function 
in the frequency plane [Q, or in terms of scattering from large radius in the Schwarzschild 
geometry g§. 

With the power-law tails in Schwarzschild spacetime well established (although a thor- 
ough understanding, especially at time-like infinity, is not complete ||), attention has been 
widened to more general situations. Tails in a Reissner- Nordstrom hole have been studied [|3| . 
Moreover, since the late time tail comes from scattering at large radius in the Schwarzschild 
case 0, it has been suggested that power-law tails would develop even when there is no 
horizon in the background [Q], implying that such tails should be present in perturbations of 
stars, or after the collapse of a massless field which does not result in black hole formation. 
In |3J , the late time behavior of scalar waves evolving in its own gravitational field or in grav- 
itational fields generated by other scalar field sources was studied numerically. Power-law 
tails have been reported for all these cases in numerical experiments, though with exponents 
different from the Schwarzschild case. These interesting results call for a systematic study 
of the late time tail in general nonvacuum and non-linear spacetimes. 



There are several interesting questions in a general analysis. Is the tail always a power 
in t? What determines the asymptotic behavior — the branch cut in the Green's function 
or the asymptotic form of the potential — and how are the two related? What determines 
the magnitude and the time dependence of the tail, and do these depend on local geometry 
and/or the presence of a horizon? 

In this Letter, we study these questions using (1) with a general V(x). Our work extends, 
and places in context, the known results for specific background geometries. 

We shall first present numerical simulations of the late time behavior for various V(x), 
showing that the decay is not necessarily an inverse power of t. Another interesting obser- 
vation is that when the parameters of the potentials are continuously varied, the behavior of 
the late time tail can change discontinuously. We then determine analytically the amplitude 
and the time-dependence of the late time tail in terms of the strength of the branch cut 
in the Green's function. The relation between the cut and the asymptotic structure of the 
potential is obtained. We establish that the local properties of the potential affect only the 
magnitude but not the time-dependence of the late time tail. 

Numerical Simulations. We study numerical evolutions of <p(x, t) given by (1) for various 
V(x). The variable x is related to, but not the same as, the circumferential radius r ||. For 
a nonsingular metric, e.g., that of a star, r £ (0, oo) maps into x £ (0, oo). For a metric 
with an event horizon at r = r Q (with g tt = at r Q Q), then r £ (r , oo) typically maps 
into x £ (— oo, oo) (the tortoise coordinate). The evolutions shown here are for the half line 
x £ (0, oo) case, with boundary conditions = at x = and outgoing waves for x — > +oo. 
The full line case [x £ (—00,00), with outgoing wave boundary conditions for \x\ — * 00] is 
basically the same. 

In the following we consider two classes of V(x): potentials which go as a centrifugal 
barrier 1(1 + l)/x 2 ( I is an integer) plus V(x), with V(x) being (i) x^~ 2 /x a ("power-law 
potentials") or (ii) (x"~ 2 /x a ) \og(x/x ) ("logarithmic potentials") when x — > +00 for some 
x a . The logarithmic potential includes the Schwarzschild metric as a special case. The 
evolution is basically independent (see below) of the initial data. The cases reported here 
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use gaussian initial data for and d<p/dt. 

Fig 1. shows log |0| versus logt at some fixed point x for several power-law potentials. 
Solid lines represent the numerical evolutions; earlier times are suppressed for clarity. After 
some quasinormal mode ringing, they approach and coincide with the analytic results (to 
be derived below) representing power law decays t _M , where /i = 21 + a except that in case 
(c) n jumps discontinuously to 21 + 2a — 2. Such jumps occur whenever a is an odd integer 
less than 21 + 3. (We assume throughout that the initial d<p/dt is not exactly zero; otherwise 
the exponents \i increase by 1, also shown below.) 

Logarithmic potentials often lead to logarithmic late-time tails. To exhibit this behavior, 
Fig. 2 shows |0|t 2/+Q versus logi for several logarithmic potentials. The numerical evolu- 
tions approach and coincide with the analytic results which represent decays in the form of 
t _M (logt) /3 , with // = 21 + a, and (3 = 1, except that (3 jumps discontinuously to (and the 
late-time tail becomes a power law) in case (c); such jumps again occur when a go through 
an odd integer less than 21 + 3. 

Table 1 summarizes these and other cases studied but not shown here. In all these 
examples, a is taken to be larger than 2 (a < 2 will be discussed elsewhere). 

Analytic treatment. We first present a heuristic picture, and then state the necessary 
modifications. Consider a wave from a source point y reaching a distant observer at x. The 
late time tail is caused by the wave first propagating to a point x' ^> x, being scattered by 
V(x'), and then returning to x, arriving at a time t ~ {x' — y) + (x' — x) ~ 2x'. Thus at late 
times oc V(x') ~ V(t/2). In particular, if V(x) ~ x~ a (logx) 13 , then one expects the late 
time tail to be ~ t~ a (logt) /3 . 

This picture requires two modifications. First, a centrifugal barrier, corresponding to 
free propagation in 3 dimensions, does not contribute to the late time tail, so that it is the 
remainder of the potential, V(x), that matters. For V(x) ~ a;~ a (logx)^, late time tail turns 
out to be t~ < - 2Z+Q! * ) (logi) /3 . The suppression by a factor t~ 21 , at least in the case a = 3, is 
known from studies of specific black hole geometries @,|3|]. Secondly, if a is an odd integer 
less than 21 + 3, the leading term in the late time tail vanishes. For (3 = 0, the next leading 
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term is expected to be t~( 2l+2a ~ 2 \ while for (3 = 1, the next leading term is £-( 2 ' +Q without 
a log t factor. 

Next we present a full analytic treatment for the half-line problem; modifications for the 
full-line case are straightforward. The evolution of <f)(x,t) described by (1) is 

t) = Jdy [G<fi(y, 0) + G<j>(y, 0)] (2) 

for t > 0, where the retarded Green's function G(x, y; t) is defined by DG = 6(t)S(x — y) 
with G = for t < and the outgoing wave boundary condition as given for above. 
The Fourier transform G satisfies 

D{u)G = [-uj 2 -d 2 x + V(x)] G{x, y- u) = 5(x - y) (3) 

and is analytic in the upper half u plane. Define auxiliary functions / and g by D(u)f(u, x) 
= D(u))g(u),x) = 0, where g satisfies \im x ^ 00 [e~ ZUJX g(u!,x)} = 1, and / satisfies f(uj,x = 0) 



= 0, f'(u,x = 0) = 1 QTU|| for the half-line problem, and \im x ^_ OQ [e luJX f(u, x)] = 1 for 
the full-line problem. In terms of / and g, and henceforth assuming y < x, G(x, y\ uj) = 
f(uj,y)g(uj,x)/W(uj), where the Wronskian W{u) = W(g,f) = g(df/dx) — f(dg/dx) is 
independent of x. Now express G in terms of G and for t > 0, distort the contour for the 
Fourier integral to a large semicircle in the lower half u plane. One therefore identifies three 
contributions to G, as follows. 

First, the integral on the large semicircle can be shown to vanish beyond a certain time 



tp(x,y) = O(x), and does not affect the late time behavior [[TT|j . 

Secondly, at the zeros of the Wronskian W(uS) at uj = ujj,j = ±1, ±2, ... on the lower 
half plane, / and g are proportional to each other, and each of them satisfies both the left 
and the right boundary conditions. Hence they are quasinormal modes (QNM's), and their 
collective contribution Gq vanishes exponentially at late times. [The case where G = Gq 



(QNM's being complete) has been discussed in detail |TT[ .] 

Finally, there may be singularities of / and g in u, which lead to the late time tail. If 
the potential V(x) has finite support, say on (0, a), then one can impose the right boundary 



condition at x = a + . Integrating through a finite distance with a non-singular equation to 
obtain g(u, x) cannot lead to any singularity in uj. It is not surprising that the same holds if 



V{x) vanishes sufficiently rapidly as x — >• +00 |IT| , |T2|| . However, if V(x) has an inverse-power 
type tail, then g(u,x) will have singularities on the negative Im uj axis (see below), in the 
form of a branch cut, as in the Schwarzschild case ||. The cut extends to to = 0, and its tip 
controls the late time behavior. For the half-line problem, f(u,x) is integrated from x = 
through a finite distance, and hence does not have any singularities in uj. For the full line 
case, / is dealt with in the same manner as g. In all cases of interest, the tail of V(x) as 
x — > —00 is either faster than any exponential, or precisely exponential. For the former, 
f(u,x) has no singularities in u, while for the latter, there will be a series of poles, but at 
a finite distance from uj = 0. In either case, the spatial asymptotics as x — > —00 has no 
bearing on the late time behavior. 

It then remains to study the spatial asymptotics as x — > +00, and the consequent 
singularities of g. First consider a power law potential with I = 0. Applying the first Born 
approximation to D(u)g = and starting with e luJX as the zeroth order solution, it is readily 
shown that 

g(u,x)~e i " x -I{u,x) , (4) 

where 

TV \ f°° j 1 s ™uj(x - x') — , iuJX , 

I(u,x)= / dx V(x )e (5) 

J x UJ 



■ . 2tuJ Xo ) a ~ 2 T(2-a) + ... , (6) 
[a - 1) L J 

The omitted terms denote a convergent power series in uj. (This form is valid only for 
nonintegral a; the integral case can be obtained by taking a limit in the final result.) The 
factor (— 2iujx ) a ~ 2 causes a cut in G on the negative Im uj axis. Thus, we have 

G (x, m t — »• 00) ~ — (7) 

v ' y ' ; #(0,0) 2 t a y ' 

The late time behavior is therefore (f) ~ t~ a generically [unless <f)(y,t = 0) = 0]. 



Notice that the time dependence is determined solely by the asymptotic form of V, 
while the magnitude, involving / and g(0,0), is sensitive to the local geometry, and hence 
to the existence or otherwise of a horizon (in the full line case). Note also that the Born 
approximation, which is strictly valid at very large x, has been used only to evaluate an 
x-independent Wronskian; therefore the results are exact at large t. 

Next consider I ^ 0. It is necessary to handle the centrifugal barrier exactly, and to 
treat only V(x) using the Born approximation. The zero order solutions are now Hankel 
functions rather than plane waves. A somewhat lengthy calculation, along the same lines as 
before, then leads to 

C (,,»;,...),- *'f"y (8) 



where 

and C(l,a) = 1 for I = 0, with F(a) = 2(2x ) a - 2 T{2l + a)/T(a), and g = 
\\m u ^Q[{iuj) l W {g , /)], which is finite |12| and reduces to #(0,0) for I = 0. The extra power 
of uj 1 in the definition of g a is responsible for the suppression of the late time tail by an extra 
factor of t~ 21 , so that in general [unless (j)(y,t = 0) vanishes], <f>(x,t) ~ f-( 2l + a ) a t j a ^ e times. 
But there is an exceptional case: when a is an odd integer less than 21 + 3, C(l, a) = and 
the late time tail vanishes in first order Born approximation; higher order approximations 
representing multiple scatterings from asymptotically far regions give the next term going 

as ~t-(2i+2a-2)_ 

Generically the late time behavior is linear in the potential (first Born approximation). 
By applying (—8/ da) on the corresponding power-law potential, we obtain from (8) that 
for logarithmic potentials: 



G(x, y; t — > oo) ~ - 



C(l,a)F{a) 

j-21+a 



9 2 o da (10) 
The leading terms are t~^ 2l+a \c\ogt + d) , except that c oc C(l, a) vanishes when a is an odd 
integer less than 21 + 3. The Schwarzschild case (a = 3) with I ^ belongs to this exception; 



as is well known 0-0], the late time behavior is a power- law with exponent — (21 + 3), and 
no logt factor. 

Analytic results for at large t can then be obtained from (7) and (8) with no adjustable 
parameters (except for case (c) in Fig. 1; see below), and are plotted as dashed lines in Figs. 
1 and 2. The agreement is perfect. The exceptional cases (c) deserve mention. For case 
(c) in Fig. 1, the leading term vanishes, and the dashed line shown represents the next 
leading term arising from multiple scattering, whose time dependence is determined, but 
whose magnitude has been left as an adjustable normalization. For case (c) in Fig. 2, the 
vanishing of the leading term implies that the asymptotic slope should be zero (i.e. no log 
t, but only a pure power, whose magnitude is determined), and this indeed agrees with the 
numerical results, with no adjustable parameters. These results are for time-like infinite. 
Results for null-infinity will be given in detail elsewhere. 

In conclusion we have achieved an analytic understanding of the late time tail in such 
systems, with asymptotic formulas agreeing perfectly with numerical results. The late time 
behavior is due to the tip of the cut in the frequency plane, which arises from scattering at 
large radius. For a potential that is a centrifugal barrier plus ~ x~ a {\ogxY(a > 2, (3 = 0, 1), 
the late time tail is generically ~ t~( 2l+a ^ (logt)^ . The possibility of a logarithmic factor in 
the leading late time behavior appears not to be widely known. Moreover, the case where a 
is an odd integer less than 21 + 3 is exceptional, and interestingly enough the most familiar 
Schwarzschild case belongs to this exception. 

We thank Richard Price for discussions. We acknowledge support from the Croucher 
Foundation. WMS is also supported by the CN Yang Visiting Fellowship and the US NSF 
(Grant No. 94-04788). 
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FIGURE CAPTIONS 

Fig. 1 log |0| versus logt, for several power-law potentials, (a) I — 0, a — 3, (b) I = 
l,a — 2.9; (c) I — l,a — 3; (d) I — l,a — 3.1. Solid lines are numerical evolutions from 
generic initial data while the dashed lines are analytical results. They are indistiguishable 
for logt > 8. To make the four sets of lines stagger, vertical shifts have been applied: (b) 
downwards by 6.0, (c) downwards by 12.0, (d) downwards by 23.0. 

Fig. 2 K\(f)\t 2l+a versus logt for several logarithmic potentials (a) I — 0, a — 3, (b) I = 1, 
a = 2.9; (c) I — 1, a — 3; (d) I — 1, a — 3.1. Solid lines are numerial results while dashed 
lines are analytic results. They are indistiguishable for logt > 9.5. For clarity, the data 
are multiplied by a constant K with (a) K = 10~ 9 (b) and (c) K = 5.6 x 10~ 10 and (d) 
K = 5.86 x 10~ 10 . 
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V(x),x — > oo 


a>2 


<f>(t),t — > oo 


x«~ 2 /x a 


odd integer < 21 + 3 


r^/i > 2/ + a 


all other real a 


£-(2Z+a) 


x*~ 2 /x a \og(x/x ) 


odd integer < 21 + 3 


£-(2Z+a) 


all other real a 


r (2i+a) logt 



Table 1 Behavior of late time tails for potentials going as 1(1 + l)/x 2 + V(x) when x — > oo. 
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